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pLn I Abstract. We derive general constraints on supersymmetric extension of axion 

(^ • models, in particular paying careful attention to the cosmological effects of saxion. 

CN , It is found that for every mass range of the saxion from keV to TeV, severe constraints 

on the energy density of the saxion are imposed. Together with constraints from axino 

f~| ■ we obtain stringent upper bounds on the reheating temperature. 
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1. Introduction 

One of the main problems of the standard model is the strong violation of CP invariance 
due to non-perturbative effects of quantum chromodynamics (QCD). In general, QCD 
effects generates the term in the lagrangian such as 

Ce = -^G" G'^''^ (1) 

where gg is the QCD gauge coupling constant, G^^^ is the field strength of the gluon, 
and G'^'^"' = e^'^'"^G" /2. Experimentally 6 must be smaller than about 10~^, but in the 
standard model there seems to be no theoretical reasons that 6 must be so small. This is 
the well-known strong CP problem. As a solution to the strong CP problem, Peccei and 
Quinn [Llj introduced anomalous global U(l) symmetry, which we denote U(1)pq. When 
this U(l)pQ is broken spontaneously, there appears a pseudo-Nambu Goldstone boson 
called axion [2j. The axion dynamically cancels the ^-term effectively, and no strong 
CP violation is observed in true vacuum. In order not to contradict with terrestrial 
experiments, astrophysical [3J and cosmological arguments |1] , the breaking scale of PQ 
symmetry Fa should lie in the range 10^° GeV ^ Fa ^ 10^^ GeV. If Fa is close to this 
upper bound, the axion is an interesting candidate for the cold dark matter. 

On the other hand, another problem in the standard model is the quadratic 
divergence of the radiative correction to the Higgs mass. In the standard model, the 
mass of the Higgs is not protected by any symmetry, and hence naturally the Higgs 
boson is expected to obtain the mass of the cut-off scale. Thus to obtain hierarchically 
small mass scale down to the weak scale requires unnatural fine-tuning. Supersymmetry 
(SUSY) [5] is the most motivated solution to this problem, since SUSY protects the mass 
of scalar fields and the weak scale becomes stable against radiative correction. From 
cosmological points of view, SUSY also provides interesting candidates for the dark 
matter. Due to the i?-parity conservation, under which standard model particles have 
the charge -|-1 and their superpartners have —1, the lightest SUSY particle (LSP) is 
stable. So if the LSP is neutralino or gravitino, they are the dark matter candidates. 

Therefore, it seems reasonable to combine these two paradigms. In fact, axion 
models are easily extended to implement SUSY. As we will see, in SUSY extensions 
of the axion models, many non-trivial cosmological consequences arise. In SUSY axion 
models, both the scalar partner of the axion, saxion, and the fermionic superpartner of 
the axion, axino, have significant effects on cosmology. But there have not been many 
studies which treat both of them in spite of their importance [El [7l [8], [9] . Moreover, these 
earlier works were based on specific models and only the restricted parameter regions 
were investigated. In this paper, we investigate all possible mass range of the saxion and 
corresponding various cosmological bounds. Our results are easily applied to any axion 
models with slight modifications, and hence provide general cosmological constraints. 
Furthermore, since the saxion and axino densities depend on the reheating temperature 
Tr after inflation, we can obtain constraints on Tr. The important result is that the 
upper bound on Tr becomes more stringent due to the late-decaying saxion and axino. 
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than that from the usual gravitino problem. 

This paper is organized as follows. In Sec. [2] we explain the dynamics of the saxion 
and its properties. In Sec. [3] we derive various cosmological constraints on SUSY axion 
models, especially upper bound on the reheating temperature. Implications for the dark 
matter in the SUSY axion models is discussed in Sec. [5l In Sec. [6] we briefly comment 
on the ultra-light gravitino scenario, where the gravitino is 0(10) eV. We conclude in 
Seed 

2. Supersymmetric axion models 

The axion is the pseudo-Nambu-Goldstone boson which appears due to the spontaneous 
breaking of PQ symmetry. The axion obtains a mass from the effect of quantum 
anomaly, but this contribution is very small. It is estimated as mao ~ 6 x 
10-'^ eV(10^2 GeV/F„). 

In SUSY extensions of axion models, the axion field forms supermultiplet, which 
contains a scalar partner (saxion) and fermionic superpartner (axino) of the axion |10j . 
The saxion mass is expected to be of the order of the gravitino mass {fn^i2)- On the 
other hand the axino mass somewhat depends on models, but generically can be as large 
as r«3/2- Although the interactions of both particles with standard model particles are 
suppressed by the PQ scale F^, they may cause significant effects on cosmology as will 
be seen in Sec. [3l 

2.1. Dynamics of the saxion 

In SUSY models, superpotentials must satisfy the holomorphy. When the real U(1)pq 
symmetry is combined with the holomorphic property of the superpotential, it is 
extended to complex U(l) symmetry, which inevitably includes scale transformation |llj . 
The invariance under the scale transformation means the existence of a fiat direction 
along which the scalar field does not feel the scalar potential. The saxion corresponds to 
such a fiat direction of the potential. But SUSY breaking effects lift the fiat direction and 
the saxion receives a mass of order 7713/2 Ijj The saxion field can develop to large field 
value during inflation, and begins to oscillate around its minimum when the Hubble 
parameter H becomes comparable to the saxion mass nis ~ ^723/2. In general such 
coherent oscillation of the saxion field has large energy density and hence its late decay 
may have significant effects on cosmology [6l El [9l [I3] . 

As an example, let us consider a model with the following superpotential: 

W = XX{^^~Fl), (2) 

where the superfields X, $ and $ have the PQ charges 0, +1, —1 respectively. The saxion 

t In gauge-mediated SUSY breaking models, the saxion receives a logarithmic potential from gauge- 
mediation effects. If the PQ scalar is stabilized by the balance between the logarithmic potential and 
the gravity- mediation effect, nis '^ '"^3/2 still holds [12] . 
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direction is easily identified as $$ = F^ with X = Olj Including SUSY breaking mass 
terms due to gravity-mediation effects, the scalar potential can be written as 

+ |A|2 {1$$ - F2|2 + |X|2 (|<|.|2 + |$|2) } , (3) 

where cx, Ci and C2 are 0(1) constants which are assumed to be positive. The potential 
minimum appears at |$| ~ |<l>| ~ Fa, but the initial amplitude of the saxion remains 
undetermined. In general, if the saxion remains light during inflation, the field value 
may naturally take the value of the order of the reduced Planck scale Mp. But the initial 
amplitude can be suppressed by introducing the Hubble-induced mass terms, which are 
induced by the saxion coupling with inflaton field through supergravity effect, given by 

Vh = H^c'j,\X\^ + c[\<!>\^ + c'^l^l^) (4) 

where c'x,c'i and C2 are 0(1) constants. If either c'l or C2 are negative, the $ or $ field 
roll away to ~ Mp during inflation. But if both coefficients are positive, the potential 
minimum during inflation is also given by |$| ~ |<i>| ~ Fa, which almost coincides with 
the low energy true minimum. Since there is a priori no reason that we expect that 
C1/C2 is exactly equals to c[/c2, these two minima are separated by 



s ~ 



C2 



Fa, (5) 



where s denotes saxion field (here we have assumed Ci > C2 and c[ > Cj). This simple 
model provides one realization of the scenario for the initial saxion amplitude Sj to 
be ~ Fa- For the case of Sj ~ Mp, the cosmological constraints become much more 
stringent than the case of Sj ~ Fa, and hence hereafter we consider only the latter case. 
For this scenario to work, the Hubble parameter during inflation Hi should be 
smaller than ~ Fa, since otherwise the large Hubble mass term during inflation takes the 
all fields to the origin. Once they are trapped at the origin, the saxion has unsuppressed 
interaction with particles in thermal bath and get a large thermal mass, which results 
in further trap of of the saxion at the origin. Thus, the oscillation epoch is significantly 
delayed |9j|§| Because delayed oscillation only makes the cosmological saxion problems 
worse, we do not consider such a case. On the other hand, the axion field has isocurvature 
perturbations with amplitude ~ Hj/^nsi) during inflation [161 IE]) which leads to a 
constraint on Hj as 

H,<2. 10' GeV .-■ [^) (I) {j^y"\ (6) 

where 6i denotes the initial misalignment angle of the axion, Q^n denotes the density 
parameter of the nonrelativistic matter and h is the present Hubble parameter in units 

I If the A-term contribution Va ~ m^^2^XFa + ^.c. is included, the X field can have the VEV 
|X| ^ m^/2- But this does not modify the following arguments. 

§ This is the case in the KSVZ (or hadronic axion) model [14,. In the DFSZ model [15] such a thermal 
mass may not arise because of the small coupling of the PQ scalar, but there arises another difficulty 
from domain wall formation. 



Cosmological implications of supersymmetric axion models 5 

of 100 km/sec/Mpc. Here we have used the observational constraint that the ratio of 
the isocurvature perturbation to the adiabatic one should be less than about 0.3 |18j . 
Thus, as long as we stick to Sj ~ Fa-, the requirement Hj < Fa is in fact valid from 
cosmological point of viewjjjj 

Although we have presented a specific model above, the dynamics of the saxion 
does not depend on axion models much, once the initial amplitude and the mass of the 
saxion are fixed. The saxion field starts oscillation at iJ ~ uis with initial amplitude Sj. 
As explained above, tUs is likely of the order of 7713/2 and the natural expectation of the 
initial amplitude is Sj ~ Fa or Sj ~ Mp. 

Now let us estimate the saxion abundance. First we consider the saxion abundance 
in the form of coherent oscillation. It is independent of the reheating temperature when 
the reheating temperature is high, i.e., F/ > rris (F/ : the decay rate of the inflaton ). 
The saxion-to-entropy ratio is fixed at the beginning of the saxion oscillation {H ~ ?7is), 
and given by 



s ) 8^°'' [mpJ ^^^ 

where Tqsc denotes the temperature at the beginning of the saxion oscillation. On the 
other hand, if F/ < m^, the ratio is fixed at the decay of inflaton H r^ Tj, and the 
saxion-entropy ratio is estimated as 



-2.1xlO-^-'GeV ;" ^tP — • (10) 

Vl05 GeVy Vl0^2 GeVy \FaJ ^ ^ 

The saxion is also produced by scatterings of particles in high-temperature plasma. 
For Tr > Td ~ lO^GeV (Fa/lO^^GeV)^, the saxions are thermalized through these 
scattering processes and the abundance is determined as [10] 



For Tr<To, this ratio is suppressed by the factor Tn/Tc- The result is 

(^Y"^' . 1.0 X 10- GeV (^^) (^^) f 15!!^y . (12) 
\sJ VI GeV/ VlO^ GeVy V Fa J ^ ^ 

Here we have assumed that thermally produced saxions become nonrelativistic before 

they decay. This assumption is valid for the parameter regions we are interested in. We 

can see that the contribution from coherent oscillation is proportional to F^ while that 

from thermal production is proportional to F'"^. Thus for small Fa thermal production 

may be dominant. Note that this expression is valid for Tp > mg. Otherwise the saxion 

II The saxion does not give rise to isocurvature fluctuation because it has large Hubble mass and its 
quantum fluctuation is suppressed during inflation. 
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log[p,/5(GeV)] 




Figure 1. Theoretical predictions for the saxion-to-entropy ratio. Thick blue lines 
represent contribution from the coherent oscillation {ps/ sp^i with Si ^ Fa and thin 
red ones represent thermal contribution [ps/sY ">. Solid, dashed and dotted lines 
correspond to Tr = 10^'' GeV,10'' GeV and 1 GeV, respectively. 



cannot be produced thermally. As a result, the total saxion abundance is sum of these 
two contributions, 

In Fig. [U we show theoretical predictions for the saxion-to-entropy ratio with 
Fa = 10^° GeV, 10^2 Qgv and 10^^ Gevfl In the figures the thick blue lines represent 
contribution from the coherent oscillation (p^/s) *■*"-' with Sj ~ Fa and thin red ones 
represent thermal contribution (p^/s)^^^^. Solid, dashed and dotted lines correspond 
to Tr = 10^° GeV,105 GeV and 1 GeV, respectively. We can see that for Fa < 
10^^ GeV {Fa < 10^° GeV), the contribution from thermal production dominates at 
Tr > IQio GeV {Tr > 10^ GeV). 



2.2. Decay of the saxion 

Since the interaction of the saxion to other particles is suppressed by the PQ scale Fa, 
it has a long lifetime and decays at cosmological time scales, which leads to several 
cosmological effects. First, let us consider the saxion decay into two axions, s -^ 2a. If 
we parametrize PQ scalar fields ^j as 



$. 



fjCXp 



.V2Fa 



(14) 



where Qi is the PQ charge of the i-th PQ field, and Fa = ^/^i(li\v^\'^, the saxion and 
axion are identified as s = Re[cr] and a = Im[cr]. The kinetic term is expanded as 

J2 I^m"^.!' ~ I 1 + ^s 1 ( ^d^ad^a + i;d^sd^s )+..., (15) 



Fa 



% Axion overclosure bound ensures O}-'' Fa < 10^^ GeV. Thus by tuning 6',;, Fa ^ 10^'' GeV is allowed. 
Late-time entropy production also makes such a large value of Fa viable, although the low reheating 
temperature Tr < 1 GeV is needed [H HOI [H] . 
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where / = ^^ q^vf / F"^. From this couphng, we can estimate the decay rate of the saxion 
into axions as 

If / ~ 1 as in many cases including the case with only one PQ scalar, this is the 
dominant decay mode of the saxion [22] . Then the lifetime is given by 

But for the model with superpotential Eq. ([2]), / can be zero at tree level due to the 
cancellation if ci = C2 in Eq. ([3]). It is crucial for cosmological arguments whether the 
dominant decay mode is into axions or not, because axions produced in the decay do 
not interact with other particles and cosmological constraints can be relaxed if this is 
the dominant decay mode. In this paper, we consider both possibilities / ~ 1 and / ~ 
and derive cosmological constraints. 

Next we consider other modes in which the saxion decays into standard model 
particles. Here we assume the saxion is lighter than SUSY particles and its decay into 
SUSY particles is kinematically forbidden. Implications of the saxion decay into SUSY 
particles are discussed in Sec. [5l For hadronic axion model, the leading contribution for 
iris ^ 1 GeV comes from the decay into two gluons. The decay rate is estimated as 

2 3 

where a^ denotes the SU(3)c gauge coupling constant. The emitted gluons produce 
hadron jets which may affect the big bang necleosynthesis (BBN) as seen in Sec. |3] 

On the other hand, the decay into two photons is always possible, which has the 
decay rate, 

r(s -. 27) ^ ^i!^!4, (19) 

where aEM denotes the U(1)em gauge coupling constant, and k is a model dependent 
constant of 0(1). These photon produced in the decay also bring about cosmological 
difficulty. 

In the DFSZ axion model, the PQ scalar has tree level coupling with the ordinary 
quarks and leptons. For m^ > 2mui{2mdi) where Ui{di) denotes the up-type (down-type) 
quark in the i-th generation [i = 1, 2, 3) the saxion decays into a fermion pair with the 
decay rate. 



■57r 



X + X ^ J V-Ta/V rni 



r(,^<(.j,)=;^(^^) ™.m (i-^l . (21) 



where x = tan/3 = {Hu)/{Hd) (through this paper, we set x = 5). Here it should be 
noticed that for rris < 1 GeV, the effective coupling of the saxion with hadrons (mesons) 
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should be used. However, there are no parameter region where the saxion decay into 
mesons has main effects on cosmology in the following discussion. Moreover, the decay 
into the muon pair (see below) gives the same order of the decay rate for m^ > 210 
MeV (Decay into mesons are kinematically forbidden for rris < 270 MeV). Therefore, 
for simplicity we neglect the effects of saxion decay into mesons for rris < 1 GeV. 

The saxion coupling with leptons is not suppressed for the DFSZ axion model, 
which gives the decay rate 

Thus we can see that the saxion decay into heavier fermions is enhanced, as long as it 
is kinematically allowed. In fact, the decay into fermions may be the dominant mode 
for some mass region even if / = 1. 

In the KSVZ model, the decay of the saxion into quarks and leptons is suppressed 
because the saxion does not directly couple with them. 

Hereafter, we consider the following four typical cases labeled as model (a)-(d). 
Model (a) denotes the KSVZ model with / = 1 and model (b) denotes the KSVZ model 
with / = 0. Model (c) denotes the DFSZ model with / = 1 and model (d) denotes the 
DFSZ model with / = 0. 

3. Cosmological constraints from saxion 

Given the decay modes of the saxion, we can derive generic constraints on the saxion 
density depending on its lifetime and mass. As we will see, for almost all the mass range 
(1 keV ^ rus '^ I TeV) the saxion density is bounded from above, although the upper 
bounds depend on the cosmological scenario such as the reheating temperature and the 
initial displacement of the saxion field. 

3.1. Effective number of neutrinos 

Relativistic particles produced by decaying particles would contribute to the additional 
radiation energy density, parametrized by the increase of the effective number of 
neutrinos, ^Ny. The definition of N,^ is given through the relation. 



Prad(T) 



7,. /T, 



I + ttAT, 



-^7 



p,(T,), (23) 



where prad denotes the total relativistic energy density, T^ and T^ denote the temperature 
of the photon and neutrino. In the standard model with three species of light neutrinos, 
N,y ^ 3.046. But if there exists another species of relativistic particle, either thermally 
or nonthermally, it contributes to the total radiation energy density parametrized by 
N^. Then the additional contribution AA^,^ is given as AN,^ = 3(pi.ad — P-y — Pu)/Pu- The 
increase of N,y speeds the Hubble expansion up and causes earlier freeze-out of the weak 
interaction, which results in '^He overproduction. The recent analyses of primordial ^He 
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abundance [231 [211 [2S] ^-^e almost consistent with Ny ~ 3. Thus we conservatively adopt 
^Ny < 1 as the BBN constraint. Note that this constraint applies to the saxion whose 
lifetime is shorter than 1 sec and whose main decay mode is s — ;> 2a |35j. 

The increase of /S.Ny changes the epoch of the matter-radiation equality and affects 
the structure formation of the universe. Thus, ANy is also constrained from cosmic 
microwave background (CMB), galaxy clustering, and Lyman-a forest. According to 
the recent analyses [271 |28l [291 [301 [31] , ^Ny ::^ 1 is not favoredjll This constraint applies 
to the saxion with lifetime Tg < 10^^ sec. 

If the saxion decays mainly into axions, AA^j, is determined from the relation, 

^ ~ 0.34(?,,(r,)-iAAr,r„ (24) 

s 

where T^ is the temperature at the decay of the saxion and g^^g counts the relativistic 
degrees of freedom. Thus the requirement /S.Ny < 1 constrains the saxion abundance as 



^ < 3.4 X 10-^ GeV ( ^^\ ( ^^ . 



(25) 



Note that almost the same constraint is applied even if the saxion decay into axions 
is suppressed, if its lifetime is longer than ~ 1 sec. The reason is as follows. Roughly 
speaking, the constraint AA^,^ < 1 means that the saxion should not dominate the 
universe before its decay. If it dominates the universe before its decay, a substantial 
amount of entropy is released by its decay. But entropy production after BBN is severely 
constrained because the baryon-to-entropy ratio should be unchanged between BBN and 
the recombination epoch, as CMB anisotropy measurements and observed abundances 
of the light elements indicates [33]. Thus in this case the constraint fl25l) is applied. On 
the other hand, entropy production before BBN is possiblelj No constraint is imposed 
in this case if the branching ratio into axions is suppressed (see the case (b) and (d) 
of Figs. [2] and [3]). The following analyses and the resulting constraints on the saxion 
abundance do not depend on whether the saxion dominates or not. 

3.2. Big-Bang nucleosynthesis 

The saxion with its lifetime > 10~^ sec may affect BBN [31]. The saxion decays into 
ordinary particles either radiatively or hadronically. If the hadronic decay occurs at 
early epoch (r^ < 10^ sec), the main effect on BBN is p ^^ n conversion caused by 
injected pions, which results in helium overproduction. At later epoch, photo- and 
hadro-dissociation processes of light elements take place efficiently. When s — > 2a is 
the dominant decay mode, the branching ratios into radiation or hadrons are small. 
Nevertheless, even a small fraction of the energy density of the saxion which goes into 
radiation or hadrons may have impacts on BBN. In particular, if hadronic decay modes 
are open, the constraint is very stringent. 

+ It is pointed out that including Lyman-a forest data raises the best- fit value of Ny [3T] , but AN^ ^ 1 
is still disfavored. If ANy is close to 1 and r^ > 1 sec, this may solve the observational discrepancy of 
Ny at BBN and structure formation [32^. 
* See Refs. [35l |36] for the case of thermal inflation driven by the saxion field trapped at the origin. 
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The constraints from photo- (hadro-) dissociation are approximately written as 

u (P£\ <\ 10-6- 10-1^ GeV for lOSec < r, < 10^ sec 
'\s) ~ 1 10-1^ GeV for lO^ec < r, < lO^^ ' l^'^J 

for radiative decay, and 

Ps\ ^ [ 10-9- 10-13 Qgv for lsec<r, < lO^sec 

sec 



^M J ^i 10-13- 10-14 GeV for lOSec < ^ < lO^^ =.- ' ^^'^^ 



for hadronic decay, where Br and B^ denote the radiative and hadronic branching 
ratios, respectively (here Br includes the hadronic decay modes). Note that if the 
injected photon energy (which is equal to the half of the saxion mass) is smaller than 
the threshold energy to destroy the light elements especially ^He, which is typically 0(10) 
MeV, the photo- and hadro-dissociation constraints are much weakened. In particular 
BBN constraints are neglected for rris ^4.5 MeV, which corresponds to the threshold 
energy for the process D + 7 — > n + p. 

3.3. Cosm,ic m,icrowave background 

The saxion with lifetime 10^ sec ^ t^ < lO^^ sec may affect the blackbody spectrum of 
CMB. Since preserving the blackbody spectrum requires the photon number-violating 
processes such as double-Compton scattering, to maintain thermal equilibrium between 
photons and electrons, photons injected in the decay distort the CMB spectrum at 
t > 10^ sec when the double-Compton scattering becomes inefficient. The distortion is 
characterized by the chemical potential yU at t < 10^ sec when the energy transfer by 
the Compton scattering is efficient, and Compton y-parameter at later epoch, which 
characterizes the deviation of the CMB spectrum from thermal distribution due to the 
inverse Compton scattering by high energy electrons. They are constrained from COBE 
FIRAS measurement as |/i| < 9 x 10"^ and y < 1.2 x 10~^ [37]. /i and y are related to 
the injected photon energy 5p^ as [3HI |39] 

^ ~ 0.714/i, (28) 

for 10^ sec < r^ < 10^ sec, and 

^ ~ 4?/, (29) 

Pi 

for 10^ sec ^ t^ < lO^^ sec. This in turn constrains the saxion energy density depending 

on its branching ratio into radiation Br, as 

1/2 



9.0 X 10-13 GeV (i^) (10^ sec < r, < 10^ sec) 

6.7 X 10-13 GeV f ^^j (10^ sec < r, < 10i3 sec) 
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3.4- Diffuse X('y)-ray background 

The two photon decay of the saxion with hfetiine longer than ~ 10^'^ sec may contribute 
to diffuse X(7)-ray background. The mass of the saxion which has such a long lifetime 
is typically smaller than 1 GeV. The photon with energy 1 keV< E.y < 1 TeV is 
transparent against the scattering with cosmic background photons and intergalactic 
medium, and hence such decay-produced photons freely propagate through the universe 
and can be observed as diffuse background photons [ID]. 

The flux of the photons from the decay of the saxion is calculated as [H] 



where 5^ denotes the branching ratio into two photons and ns{z) is the number density 
of the saxion at the redshift z. E' is the energy of the photon at the instant of 
production and E is the present redshifted energy, the relation between them is given 
by £" = (1 + z)E. Under the assumption of the fiat universe (^a + Q^ = 1), this 
expression can be integrated yielding 






X exp 



lJv^^(fl)-l)(v^+l) 



('32) 
_3i7or.v^"^(v^(?(f|) + l)(v^-l) ^ ^ 

where n^o denotes the present number density of the saxion, Hq denotes the present 
Hubble constant, and 

g{x)= [^A + fi^x^]-'/'. (33) 

On the other hand, the observed photon flux in the range 1 keV ^ E < 100 GeV 
is roughly given as 

r 8{^)-'-' (0.2 keV < E < 25 keV) 

F,oUE) r^l 57x 10-^ ( Jv)"''' (25 keV < i^ < 4 MeV) (34) 

[ 17 X 10-6 (hk™)"'' (4 MeV < i^ < 120 GeV), 
in the unit of cm^^ sec^^ sr"^, from the observations of ASCA [42], HEAOl [33], 
COMPTEL [44J, and EGRET |35]. Thus, from the requirement F^ < F^^hs, the tight 
constraint on the saxion density is derived. Note that even for the saxion lifetime 
longer than the age of the universe, the small fraction of the decayed saxion at t < to 
contributes to the diffuse background and its abundance is limited. 

To estimate this constraint, let us consider the case with r^ > tg- In this case, the 
X(7)-ray spectrum of photons from the saxion decays has the maximum at -Emax = "^s/2, 
where the flux is simplified as F^(ii^max) = B^Usq/ iflnTsHo) . Then from the condition 
-^7(-Emax) < F^ohs{Ems.y), wc obtaiu a constraint, 

B, (-j < 2vr^— F,.,, 1^- 



2.4/i X 10"'" GeV 



mo 



1 MeVV 
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X ( ^' \ ( -^7°bs(»^s/2) \ 

V 1018 sec/ V 10-2 cm-2sec-V ' 
where sq denotes the present entropy density. For the case with r^ < to, the photon 
energy which gives the flux maximum deviates from E = ms/2 due to the Hubble 
expansion. For r^ <C to, it is given by -Emax = (ms/2)(3i/oTsV^m/2)2/'^. This leads to 
the constraint 

fPs\ 47rm, 

-D7 I j ~ ^ -^7obs l-f^maxj 

V S / <J Sg 

~ 4.8 X 10-19 GeV (^^) f ^^^^}dE^pL\ . (36) 



cJ.J. Reionization 



If the saxion decays after recombination era and the injected photon energy is relatively 
small (ms ^ 0(1) keV-O(l) MeV and 10^^ sec < r^), redshifted photons may leave 
the transparency window until the present epoch [IQ]. Then, emitted photons interact 
with and ionize the intergalactic medium (IGM), and they contribute as an additional 
source of the reionization. If this contribution is too large, the optical depth to the 
last scattering surface is too large to be consistent with the WMAP data [33j. Here 
we apply the results from Refs. |10l US] , simply assuming that if the decay-produced 
photon leaves the transparency window, one-third of the photon energy is converted 
to the ionization of the IGM (the remaining goes to the excitation and heating of the 
IGM). According to Refs. [lOl HH], this is a good approximation when the decay occurs 
before the reionization due to astrophysical objects takes place and most of hydrogen 
atoms exist in the form of neutral state. (The Gunn-Peterson test indicates that the 
reionization occurred at z ~ 6 [47J.) The constraint on the saxion density can be written 
as 

ft (^) < (^) . (37) 

V S / \ S / bound 

where (Ps/s)^,^^^^ can be read off from Fig. 2 of Ref. |46]. For example, for r^ > to, it is 
given by 

(^) - 4.3 X 10-1^ GeV (^^) (^ , (38) 

V s /bound V1018 sec/ V 0.0227 ^ ^ 

where Qb denotes the density parameter of the baryonic matter. This constraint is 
complementary to the diffuse X(7)-ray limit. 

3.6. Present matter density limit 

For the saxion with its lifetime r^ > to, its energy density contributes to the dark matter 
of the universe, and hence the saxion density should be less than the observed matter 
density, flsh"^ ^ ^mh"^- In terms of the saxion-to-entropy ratio, this is written as 

^ < 4.7 X 10-1° GeV (^] . (39) 
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3. 1. LSP overproduction 

If the saxion mass is larger than about 1 TeV, the saxion can decay into SUSY particles. 
Here we suppose that the LSP is the hghtest neutrahno. The decay into SUSY particles 
were investigated in detail in Ref. [18] and it was found that decay into gauginos has 
roughly the same branching ratio as that into gauge bosons. Thus we should be careful 
about LSP overproduction from the saxion decay. The resultant abundance of the 
LSP depends on T^, and for T, > mLSp/20 LSPs produced from the saxion decay are 
thermalized and have the same abundance as that expected in the standard thermal relic 
scenario of the LSP dark matter. In this case no upper bound on the saxion abundance 
is imposed. On the other hand, if Tg < mLSp/20, the abundance of the LSP is given by 



Plsp 



^^--LSP^ + ^^i^^L for nLSPm)M < HiTg), 
nig s s 

87r^g^{Ts) {av)TsMp 

where B^ denotes the branching ratio of the saxion into SUSY particles, {av) denotes 
the thermally averaged annihilation cross section of the LSP, and Pl's™^' denotes the 
contribution from thermal relic LSPs taking account of the dilution from the saxion 
decay. The LSP number density immediately after the saxion decay nLSp(^s) is defined 
as TiLSP = '2BsPs{Ts)/ms. For deriving the constraint, we ignored the contribution to 
the LSP production from thermal scattering processes. Moreover, the second line of 
Eq. (l40l) always results in overproduction of LSPs with the annihilation cross section for 
ordinary neutralino dark matter. Thus for deriving the constraint, we consider only the 
first term of the first line of the right hand side of Eq. PUI) . The bound can be written 
in the form 

^ < 2.4 X 10-° GeV (^ (^^) , (41) 

for nis ^ 1 TeV and Tg < ?71lsp/20. This constraint can be relaxed if the annihilation 
cross section of the LSP is significantly large. We will revisit this issue in Sec. 
Hereafter we set mLSP = 500 GeV as a reference value. 

Including all of these constraints, we can derive general upper bounds on the saxion- 
to-entropy ratio as a function of the saxion mass nig for models (a)-(d). In Fig. [2llH we 
show the results with Fa = 10^°, 10^^ and 10^^ GeV, respectively. In each panel, the 
orange line represents the bound from AA^,^ < 1, the thick-solid brown line represents 
the bound from BBN, the thick-dotted purple line represents the bound from CMB, 
the thick-dot-dashed green line represents the bound from diffuse X(7)-ray background, 
the thin-dot-dashed blue line represents the bound from reionization, the thin-dashed 
red line represents limit from the present matter density, and the thick-dashed gray line 
represents LSP overproduction limit from the saxion decay. We also show the theoretical 
prediction for the saxion energy density in the figures for Tf> = 10^*^ GeV and 1 GeV by 
thin-dotted black lines. 
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Figure 2. Various cosmological constraints on the saxion abundance for Fa = 
10^" GeV. Thin dotted black lines represent theoretical prediction p^/s = (ps/s)''*^'' + 
{pjs)^^^^ for Tr ^ IQi" GeV (upper) and Tjx = 1 GeV (lower) with Si = Fa. Four 
panels correspond to different models. Model (a) : KSVZ with f ~ I, model (b) : 
KSVZ with / = 0, model (c) : DFSZ with / = 1, model (d) : DFSZ with / = 0. 



4. Constraints from axino 

So far, we have ignored the cosmological effects of axino, the ferniionic superpartner of 
axion. The mass of axino is model dependent, but it can be as heavy as the gravitino 
mass, 7713/2 lini [22l [50]. For example, in the model of Eq. ([2]), the axino mass is 
estimated as rria ~ |AX|. As noted earlier, taking into account the A-term potential like 
Va ~ m3i2\XF'^ + h.c, X can have the VEV of the order of 77^3/2 • Thus in this model 
the axino mass is naturally expected to be 7713/2 • Hereafter for simplicity we assume 



ma, ~ 



"^3/2- 



Axinos and gravitinos are produced through scatterings of particles in thermal 
bath. First, we assume either of them is the LSP, and hence their thermally produced 
abundance must not exceed the present abundance of the dark matter [Sljl^ The 

tt As long as ttIj ^ m.3 /2 ~ m^, which of them is lighter is not relevant. But if ma <gC TO3/2 and m-^M ^ 
1 GeV, thermally produced gravitinos can decay into axinos within the present age of the universe and 
the constraint is relaxed [52]. 
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Figure 3. Same as Fig. [21 but for Fa = lO^^ Q^y 



abundance of thermally produced axinos is calculated as [53] 



Pa 



on Tn-7 6 n ^r f ^'^ WiU"LreV\ / 1r 



(42) 



where (7^ is the QCD gauge coupling constant. Thus for large ma(~ m^), the constraint 
on the reheating temperature becomes stringent. Note that axino thermal production 
for Tr < ITeV is negligible, because SUSY particles are not produced efficiently for 
such a low reheating temperature. 

Next, if the axino mass is larger than about 1 TeV, the axino can decay into SUSY 
particles. This LSP abundance depends on the temperature at the decay of the axino, 
Ta similar to the case of the saxion decay. If T^ > mLSp/20, LSPs produced from the 
axino decay are thermalized and the standard thermal relic scenario of the LSP dark 
matter is maintained. If Ta < mLSp/20, the LSP abundance is determined by Eq. (HOll 
after replacing Tg and ps/s with T^ and Pa/s. 

ff Here we assurae there is no entropy production after the reheating ends. If the saxion dominates 
the universe and decays before BBN, the axino and gravitino abundance can be reduced. 
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Figure 4. Same as Fig. H but for Fa = 10^'' GeV. 



On the other hand, the thermally produced gravitinos have the abundance as 

■lGeV\ / ms \2 / Tr 



TP 
Pzl2 



6.3 X 10"^' GeV 



m-i/2 



-\ 



ITeV/ V 106 GeV 



(43) 



for m3/2 <^ mg where nig denotes the mass of the gluino (here the logarithmic 
dependence on Tr is omitted). Contrary to the axino, the constraint becomes severer 
when wi3/2(~ Tns) becomes smaller. Note that for m^i2 ^ 1 keV, gravitinos get 
thermalized and their abundance becomes independent of the reheating temperature. 
But for 16 eV < m^/2 ^ 1 keV they contribute to the dark matter density as a 
hot component due to their long free-streaming length, and such a contribution is 
constrained from cosmological observations in particular Lyman-a forest data [55] . Thus 
the gravitino mass in this region is strongly disfavored. (The case of ultra-light gravitino 
m^j2 < 16 eV will be mentioned later.) In addition to the constraint from the present 
matter density, there may be another constraint coming from the late-decay of SUSY 
particles into gravitinos or axinos, which may affect BBN. But the constraint is quite 
model dependent, and hence we do not consider it here. 

Including those constraints from the gravitino and axino, we derive the upper 
bounds on the reheating temperature for each saxion mass and show them in Figs. [MZl 
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Figure 5. Upper bounds on the reheating temperature Tfj for each model with 
Fa = 10^*^ GeV. The initial amplitude of the saxion is assumed to be Si ^ Fa- The 
thin-short-dashed light blue line represents the bound from axino thermal production, 
and thin-dotted black line represents the bound from gravitino thermal production. 
The shaded region contradicts with the lowest possible reheating temperature [56j . 
The other lines are the same as Fig. [2j Four panels correspond to different models. 
Model (a) : KSVZ with / = 1, model (b) : KSVZ with / = 0, model (c) : DFSZ with 
/ ^ 1, model (d) : DFSZ with / = 0. 



In these figures, we have assumed that the initial aniphtude of the saxion is given by 
Si ^ Fa- The thin-short-dashed light blue line represents the bound from axino thermal 
production, and thin-dotted black line represents the bound from gravitino thermal 
production. The other lines are the same as Fig. [2l Because the saxion-to-entropy 
ratio (ps/s)*^^) is proportional to (Trs^), the constraints on Tr scale as s~^. However, it 
should be noted that although the axino constraint is stringent for relatively large m^ 
as can be seen from these figures, it may be significantly relaxed if the axino is much 
lighter than the gravitino. 
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5. Dark matter candidates 

We have seen that the very stringent bound on the reheating temperature is imposed for 
wide range of the saxion mass. It is typically stronger than the usual upper bound from 
the gravitino overproduction. It has some implications to the dark matter candidates. 
For example, it invalidates the gravitino dark matter for wide parameter regions. In 
this section we summarize the dark matter candidate in SUSY axion models. 



5.1. TTis < 1 TeV 

First consider the case where the gravitino (axino) is the LSP. In our model, we assume 
the axino has the mass comparable to the gravitino, and hence the axino (gravitino) is 
the NLSP. Which of them is the lighter is not important because both have the similar 
properties. The saxion also has the mass comparable to the gravitino. The stability 
of the saxion is not ensured by the i?-parity, but since its decay rate is suppressed by 
the PQ scale Fa, the saxion lifetime can exceed the present age of the universe. If this 
is the case, the saxion can be the dark matter. In addition, if 9]''^Fa ~ 10^^ GeV, the 
axion can also play a roll of the dark matter. Therefore, we have four candidates for the 
dark matter, i.e. gravitino, axino, axion and saxion. However, the saxion dark matter 
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Figure 7. Same as Fig.[5l but for Fa = IQ^^ GeV. 



is possible only when Fa ~ 10^^ GeV and 1 keV< rus ^10 keV, as can be seen from 
Fig. [71 and hence is less attractive candidate. For 10^° GeV < -Fa ^ 10^^ GeV, the 
axino dark matter is allowed for wide parameter regions, but the gravitino dark matter 
is excluded except for rus < 100 keV, as can be seen from Fig. [5] and Fig. [61 On the 
other hand, for Fa > 10^^ GeV, the axino or gravitino dark matter is almost impossible 
and the axion is the most viable dark matter candidate, although other cosmological 
constraints are severe (Fig. [7]). 

To summarize, for irts ^ 1 TeV, the axino is a good dark matter candidate for 
Fa < 10^^ GeV, and the axion may be dark matter candidate for Fa > 10^^ GeV. 



5.2. nis > 1 TeV 

For larger rrisi— iTL^/2), the saxion decay mode into SUSY particles opens. Then the 
LSP is assumed to be the lightest neutralino. As discussed in Sec. 13.71 the abundance is 
given by Eq. ( HOl) if the decay occurs after the freeze-out of the LSP. For the neutralino 
which has small annihilation cross section, very low reheating temperature is needed to 
obtain the correct abundance of the dark matter as can be seen from Figs. [MZl However, 
for the neutralino with larger annihilation cross section such as wino- or higgsino-like 
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LSP, the abundance may be significantly reduced and they can become the dark matter 
independent of the reheating temperature jUj. This is easily seen from Eq. (HOj) . which 
describes that the resulting abundance is bounded from above. Thus for the saxion mass 
larger than 1 TeV, these non-thermally produced neutralinos may be dark matter. On 
the other hand, if T, is larger than the freeze-out temperature of the LSP Tf ~ m-LSp/SO, 
the standard thermal relic scenario holds and the saxion has no impact on cosmology 
unless it dominates the universe before the decay. If the saxion once dominates the 
universe, it must not decay into axions mainly, since otherwise AA^^, ^ 1 holds at BBN 
and contradicts with the observation. 

Therefore, for irts > 1 TeV, the lightest neutralino produced either thermally or 
non-thermally is likely dark matter depending on irts and Fa- Of course, the axion is 
also a good candidate for the dark matter for Fa > 10^^ GeV. 

6. Ultra-light gravitino scenario 

We have seen that for almost all the mass range of the gravitino, the reheating 
temperature is severely constrained. However, for an ultra-light gravitino with mass of 
the order of 1-10 eV, gravitinos are thermalized and their abundance is sufficiently lower 
than the dark matter abundance, and hence the reheating temperature is not constrained 
from gravitino overproduction. Thus thermal leptogenesis using right-handed neutrino 
[57] . which requires Tr ^ 10^ GeV, may be possible. 

In SUSY axion model, the saxion oscillation also contributes to the present dark 
matter density as given in Eqs. (IHl) and ( TTOj) . We should ensure that this contribution 
does not exceeds the present matter density of the universe. From Eq. (jHj), pa/ s is 
bounded as 

Thus for Fa ^ 10^^ GeV, the saxion abundance from coherent oscillation is smaller 
than that of the dark matter. However, we should be aware that thermally produced 
saxion and axino abundances are comparable to that of the gravitino for Tr > 
10^GeV(F(j/10^^GeV)^ if m^ = m^ = m3/2, since both are thermalized in the early 
universe. Thus constraint on the gravitino mass may become more stringent by a factor 
of two, if thermal leptogenesis is assumed to work. In Fig. [HI we show allowed parameter 
regions on F^-Tr plane for iris = 10 eV (indicated by (a)) and for nis = 1 keV (indicated 
by (b)) with an assumption Si = Fa- It can be seen that for rris ^10 eV, thermal 
leptogenesis is still possible as long as the PQ scale satisfies Fa '^ 5 x 10^^ GeV. 

7. Conclusions and discussion 

We have investigated the cosmological constraints on supersymmetric axion models 
which are motivated from particle physics point of view. It is found that the presence 
of the saxion and the axino makes it rather difficult to construct a viable cosmological 
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scenario which is free from any contradiction with observations. In particular, in almost 
all range of the gravitino mass (which is assumed to be the same order as the saxion 
mass), the strict upper bound on the reheating temperature is imposed, and it is more 
stringent than the bound from the usual gravitino problem [Figs. [5J1Z]- The constraint 
depends on whether the main decay mode of the saxion is into axions or not. It should 
be noted that although the axino constraint is stringent for relatively large m^ as can be 
seen from our results, it may be significantly relaxed if the axino is much lighter than 
the gravitino. The axion is a good candidate for the cold dark matter, although the 
axino or gravitino dark matter is also viable for some parameter regions. 

The obtained stringent bound on the reheating temperature has some implications 
on the baryogenesis scenarios. As is well known, the standard thermal leptogenesis 
scenario using right-handed neutrino [57] is incompatible with the gravitino problem 
except for 7713/2 ^ 10 eV or 7/13/2 ^ 10 TeV. Although the gravitino mass around 
?7i3/2 ~ 10 GeV may also be compatible with thermal leptogenesis, in this region the 
NLSP decay into gravitino may cause another difficulty The presence of saxion makes 
this situation worse, and hence the standard thermal leptogenesis does not seem to 
work in supersymmetric axion models. The Affieck-Dine baryogenesis scenario |58] can 
work well even for such a low-reheating temperature [59], except for the case of gauge- 
mediated SUSY breaking models for small 7713/2, where the Affleck-Dine mechanism may 
suffer from Q-ball formation [60]. On the other hand, for an ultra-light gravitino mass 
"^3/2 ^ 10 eV, thermal leptogenesis is still possible. Here it should be noticed that these 
constraints from the saxion strongly depend on the initial amplitude of the saxion Si and 
we have assumed Sj is roughly given by the PQ scale Fa- A concrete example which gives 
such initial amplitude is given in Sec. 12.11 Perhaps this is the smallest value expected 
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from naturalness, and hence our bounds presented in this paper should be regarded as 
conservative, which means our constraints cannot be relaxed without significant changes 
in the cosmological scenario, such as additional late-time entropy production. 

Finally we comment on the detectable signature of the SUSY axion models. If 
the saxion mass is around a few MeV, its decay into electron-positron pair and their 
annihilation may be observed as the 511 keV line from the Galactic Center. Such line 
gamma photons are actually observed [6TJ. However, in order to explain the observed 
fiux, huge entropy production that dilute the saxion density is needed [62]. Also it may 
be possible that the nature of the axion sector is determined from collider experiments 
if the axino is the LSP and charged particle such as stau is the NLSP [63] . 
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